// Copyright (C) 2009,2010,2011,2012  Marco Restelli
//
// This file is part of:
//   LDGH -- Local Hybridizable Discontinuous Galerkin toolkit
//
// LDGH is free software: you can redistribute it and/or modify it
// under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// LDGH is distributed in the hope that it will be useful, but WITHOUT
// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
// or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public
// License for more details.
//
// You should have received a copy of the GNU General Public License
// along with LDGH. If not, see <http://www.gnu.org/licenses/>.
//
// author: Marco Restelli                   <marco.restelli@gmail.com>

// Set the mesh size
h = 1000.0;

// Build the mesh
Point(1) = {    0,    0,0,h};
Point(2) = {25600,    0,0,h};
Point(3) = {25600,25600,0,h};
Point(4) = {    0,25600,0,h};
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,1};
Line Loop(5) = {1,2,3,4};
Plane Surface(6) = {5};
Extrude {0,0,6400} {
  Surface{6};
}

// This already defines the global volume
Surface Loop(29) = {27, 6, 15, 19, 23, 28};

